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Abstract 

We evaluate the correlation function of the spectral staircase and use it to evaluate the mesoscopic 
particle number fluctuations in integrable systems. 
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I. INTRODUCTION 

In the preceding paper [l], we addressed the question of mesoscopic fluctuations vis-a-vis 
thermal fluctuations in classically integrable and chaotic systems. Mesoscopic fluctuations 
can be expressed in terms of the semiclassical expression for the level density correlation 
function and thermal fluctuations in terms of the fluctuations of thermal occupancy. We 
illustrated our formalism for both classically chaotic and integrable systems by calculating 
the particle number fluctuations and fluctuations of the speciflc heat in a two-dimensional, 
non-interacting electron gas. The results are particularly interesting in the classically inte- 



grable case where the ensemble of hard-wa. 
ratios was considered as a model system 



1 rectangles with equal areas but varying aspect 
The key feature of such systems is the near 
absence of correlations between energy levels, with mean level spacing A, on scales less then 
~ \/eA in proximity of level energy e (Fermi energy e^? in this case), and strong corre- 
lations on larger scales. Level correlations were descried in terms of the correlation function 
of level density using a simplifled ansatz, which neglects their long-range oscillatory behav- 
ior but is valid at flnite temperature. As a result, we found that the speciflc heat 
fluctuations grow linearly with temperature oc T/A for T <^ Em and fall off exponentially 
for T ^ Efn- For the particle number fluctuations at a flxed chemical potential (which is 



exponentially close to Fermi energy), the zero-temperature contribution oc ^Jsp/ ^ comes 
from the variations of the number of levels over the Fermi sea. For T <^ Em-, the leading 
temperature-dependent term in mesoscopic fluctuations is negative, oc —T/A, and cancels 
out the fluctuations due to fluctuations of thermal occupancy, oc T/A. As a result, the over- 
all temperature-dependent part of the fluctuations is quadratic, oc T^/ (i^^A). For T ^ E^, 
on the other hand, the temperature-dependent part of the fluctuations is dominated by the 
fluctuations due to thermal occupancy fluctuations, oc T/A. 

In this work, we reexamine the mesoscopic particle number fluctuations. Towards this 
end, we derive a general semiclassical expression for the correlation function of the spectral 
staircase and show that it can be expressed concisely in terms of the level number vari- 
ance. It is equivalent to the semiclassical expression for the level density correlation function 
[l| but lends itself more readily to thermal averaging in calculations of the fluctuations of 
thermodynamic quantities. We use this expression to evaluate the particle number fluctu- 
ations in a two-dimensional, non-interacting electron gas in a rectangular box and flnd the 
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results consistent with those in which utihzed a simphfied ansatz for the level density 
correlation function. 



II. CORRELATION FUNCTION OF SPECTRAL STAIRCASE 

In the quasi-continuous limit, e ^ A, the spectral staircase is defined as 

Af{e)= f dep{e) (1) 



where p {e) is the level density. It is well known that rescaling the energy variable e {e) 
allows for universal Representation of the staircase A/" (e) = ep that eliminates the par- 
ticular shape of AT [5] and reduces the spectrum into that of an infinite two-dimensional 
system with constant p. Consequently, in what follows we limit our consideration to 
the thus reduced spectrum. Further below, we make an approximation that its prop- 
erties are applicable to the size-quantized 2D systems, such as a hard-wall wells, where 



1 + ^^A/ej ; 6]. neglecting corrections, we have 



A 1 

and 77 {s) = e/ A, where A is the system area. 

The correlation function of spectral staircase can be written using semiclassical expression 



for the staircase 



as 
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m^^mw) = E I ) = w (sf - E ^ (^) 



(3) 



where 6Af = Af—J\f, hw = 62 — Si <^ e = (£2 + ^i) /2, 2N is the dimension of the phase 
space, Aj (e) and Tj (e) are the amplitudes and periods of the periodic orbits, and 

^=^E^ (4) 
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Over-bars are used above to denote ensemble averaging. Using notations of 
be rewritten as 

1 



eq. can 



m{e^) 5M{e2) = - [TT {e) - S {e, \hw\)) (5) 

where. S (e, l/zo'l) is the level number variance on the interval \fvjj\ and [e) is the satu- 
ration level number variance, averaged over the oscillations. 



A quick check confirms that 



dsi de2 



5U (e,) 5U (e,) = ^ ^? (^) = (^i) (^2) (6) 

where Sp = p — p. In what follows, we set /i = 1. Generally, the level density correlation 
function can be written as 



6p{e^)6p{e2) = -6{co)-IC{e,uj) (7) 

where 5-function is universally present and indicates absence of level correlations and /C {e, u) 
describes level repulsion, is system-specific and becomes important for u greater than the 
scale for onset of level rigidity. Per ([5]), the spectral staircase correlation function reduces 
to the properties of the level number variance. For small \uj\, the universal result for S is 81] 

S(^,|^|) = ^ (8) 

and is equivalent to the (5-function term in the correlation function of level density ([7]), as is 
also seen from ([5]) and IQ. For integrable systems, for large > Em, S (e, oscillates 
persistently around S (e), and for a rectangular box is shown in Fig. 4 of Per 
([5]), this proves corresponding oscillations around zero of the correlation function of spectral 
staircase, which, in turn, points to the persistence of correlations for |ci;| ^ E^,- 

III. THERMODYNAMIC MESOSCOPIC FLUCTUATIONS. 

As was shown in pL,], the fluctuations of a thermodynamic quantity G that can be ex- 
pressed as an integral over the spectrum 

G = J dep{e)f{e)g{e),G = j de-p {e)J {e) g {e) (9) 

can be separated into mesoscopic fluctuations, due to sample-to-sample variations of identi- 
cally prepared systems, and fluctuations due to the fluctuations of thermal occupancies: 



5G^ = 6GI + 6G} (10) 




6Gl= deide26p{ei)6p{e2)f{ei)f{e2)g{ei)g{62) (11) 



dG\^ fdN\ ^ {dG\^T^ 
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where / (e) is the thermal occupancy, /i is the chemical potential, C is the specific heat, is 
the particle number in the system, and 6G = G — G. In what follows, we consider only Fermi 
thermal occupancy since we are interested in properties of a non-interacting electron gas in 
a 2D hard-wall potential. In 2D, /i and Sp are exponentially close, ^-Sp = — Texp {—ep/T), 
and the fiuctuations of thermal factors can be easily evaluated using jsj 

— ^ A, C ~ — — 13) 

and, in particular, 

T 



5N} ^ -, 6Cj ~ C (14) 

for the fiuctuations of the number of particles and specific heat.[l| 

Further, integration of (fTTj) by parts, with the use of (j6]), yields an alternative form of 



mesoscopic fiuctuations 5G'j, 

6Gl = 1 1 d6id626X (£i) 6X {82) ^•^[^^^ ^•^[^'^ g (ei) g (^2) (15) 

Notice, that (fTTj) and (fTSll represent a combined ensemble and thermal averaging. Further 
manipulations with Fermi factors yield 

a7(.097M_^-.„(l^ -A (16) 



dei de2 \ T ' 2T^ 

H{u,v) = [2 {ucschu + vcschv)]^ (17) 



u)du = ^ (18) 



where Ti. has the following properties: 

H {u,v) dv = - {—1 + ucothu) cosh"^ u, / h 

-00 2 J — CO 

Substituting ([3]), ([HD and ([I7D into ^ yields 

SGl = l I dedJ- (E- (.) - E (., I^D) " |) ^ + |) (19) 

rferfcui (S- (.) - S (., IcD) T-^n (i^, ^) ^ - I) + I) (20) 
In particular, substituting g = 1, we find 

SN', = I dsYT is) h (i^) I dedu (S (., |..|)) T^^U (i^, (21) 

for mesoscopic particle number fiuctuations. 



In terms of the zero-temperature fluctuations 



the latter can be rewritten as 




deduj{E{e, \uj\))T-^n 



e — Bp ijJ 



T ' 2T 



+ 



(22) 



(23) 



Given that structure of the function 7i, S {e, can be approximated by S [sp, and 
we flnd, neglecting the terms of order T^/e|. 

/oo 
duiJ:iep,2T\u\))hiu) (24) 
-oo 

We proceed to apply this expression for a particular form of S in a rectangular box.| 



IV. FLUCTUATIONS IN A 2D ELECTRON GAS IN A RECTANGULAR BOX 



Rectangular box is a particular case of a hard-wall potential, whose main distinction is 
that it represents a "generic" classically integrable system, that is, with no degeneracies 
that may be caused by extra symmetries. The ensemble is understood as a collection of 

1/2 

rectangular boxes of equal area A = L1L2 and varying aspect ratio aasp = -^^2/-^i-[2| this 
case, we have 



E(£,E) 



00 00 



sm 



TT^A 



Mi=0 A/2=0 



M" 



EJin/eA) (^Mfalil + M^aasT 



1/2 



3/2 



e 

vr^A 



00 00 



A/i=0 Af2=0 



sm 



E^i-K/eA) {Ml + M 



(M2 + M| 



2 A 3/2 



(25) 



(26) 



TT^A 



CXD 00 



2(5- 



M 



EE — 



3/2 



TT^A 



00 cx> 

E E 

A/i=0 A/2 = 



25- 



M 



(Mf + M| 



2\3/2 



(27) 
(28) 



Ml = 0, M2 = 
6m = 1/4 Mi,2 = 0,M2,i ^ 

1 Ml 7^ 0, M2 ^ 

where for simplicity we assumed eta™ ~ 1, as is customary in analytical evaluations of the 
sums and numerical calculations. {4], p] 







Numerical evaluation of the sum in (127|) is straightforward and we find 

(SN^) =?!lM~l^,/££ (29) 



Also, in the limit T <^ i^^, where ([8]) can be used,[l| we find 

- i^K)o - r du^-^h (u) = {SN^X - I (30) 

in complete agreement with 

Insofar as the integral in ^IM is concerned, we could not find its closed form. However, 
a very good approximation for (ITSl) can be obtained by replacing 

^(w)^— — TT — r, (31) 
vr^ smh [ua) 

For instance, a k, 2 describes its asymptotic behavior at -u — > oo, while a ~ 3/2 is suitable 
for the origin. Using (!3T]) . the integral in fl2^ can be easily evaluated and we find 



,2 



oo oo 



To proceed further with the evaluation of the sum, we must now consider two limits: T <^ 
and T ^ E'm- In the first limit, we can introduce continuous variables x,y cc {T / E„^ Mi^2 
and approximate summation with integration, which, in polar coordinates, yields a con- 
vergent integral rfp (tanhp/p)^. Consequently, we find ^A*"^ ~ i^-^p) — const {T/ A), 
as was already obtained above (l30l) .(lol| In the opposite limit, tanh 1 and we find 
~ i^-^p) ~ const {Em/ A), again, in complete agreement with 
It must be pointed out that the results in [l| were obtained using a simplified ansatz for 
the level density correlation function 

S (uj) sm{27[uj/E„ 



Sp (e) 5p{e + uj)^^ \ ' (33) 

Correspondingly, it is clear that in the finite-temperature limit the oscillations of S (e, E) 
as a function of E, per ( |25l) . can be neglected and a simplified ansatz 

T.(,M)='4d04-\ (34) 



A V E, 



■'m 



D(x)^l-2 ""'^""^ - vrxSi(2.x) 



where Si is the sine integral, can be used, which yields results identical to those in 



V. SUMMARY 



We derived the correlation function of the spectral staircase in terms of the level number 
variance, eq. We also showed that the oscillatory behavior of the level number variance 
can be neglected at finite temperatures and the simplified ansatz (IMI) can be used for 
evaluation of mesoscopic fluctuations of thermodynamic quantities, such as particle number 
fluctuation. 
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